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We show that new massless Dirac fermions are generated when a slowly varying periodic po- 
tential is applied to graphene. These quasiparticles, generated near the supercell Brillouin zone 
boundaries with anisotropic group velocity, are different from the original massless Dirac fermions. 
The quasiparticle wavevector (measured from the new Dirac point), the generalized pseudospin vec- 
tor, and the group velocity are not colhnear. We further show that with an appropriate periodic 
potential of triangular symmetry, there exists an energy window over which the only available states 
are these quasiparticles, thus, providing a good system to probe experimentally the new massless 
Dirac fermions. The required parameters of external potentials are within the realm of laboratory 
conditions. 



Semiconducting and metallic superlattice structures 
are now routinely used in manipulating the electronic 
structure of materials [ij. These superlattices have ad- 
ditional electronic band gaps at the supercell Brillouin 
zone (SBZ) boundary, which often give rise to interest- 
ing phenomena. 

Since the successful isolation of graphene 0, S 0, 0] , 
numerous studies have been performed on this novel ma- 
terial 0. In particular, there have been a number of 
interesting theoretical predictions on graphene superlat- 
tices (defined to be graphene under an external periodic 
potential or graphene with periodic defects). For exam- 
ple, for an one-dimensional (ID) or a two-dimensional 
(2D) rectangular graphene superlattice, the group ve- 
locity of the low-energy charge carriers is renormal- 
ized anisotropically [7]; a corrugated graphene sheet is 
expected to show charge inhomogeneity and localized 
states [8] ; and arrays of anti-dots (missing carbon atoms) 
of specific design could induce band gaps [Qi] or mag- 
netism [io|. 

Graphene superlattices are not only of theoretical in- 
terest, but have also been experimentally realized. Super- 
lattice patterns with periodicity as small as 5 nm have 
been imprinted on graphene through electron-beam in- 
duced deposition of adatoms (llj . Also, triangular pat- 
terns with ^10 nm lattice period have been observed for 
graphene on metal surfaces [11, [l^, 14 1 . Using periodi- 
cally patterned gates is another possible route to make 
graphene superlattices. 

In this paper, we show that when a periodic potential is 
applied to graphene, a new generation of massless Dirac 
fermions is formed at the SBZ boundaries. The elec- 
tronic wavevector (measured from the new Dirac point), 
the group velocity and a generalized pseudospin vector, 
defined below, of the newly generated massless Dirac 
fermions are not coUinear anymore. In ID or 2D rect- 
angular graphene superlattices, the features of these new 



massless Dirac fermions are obscured by other states ex- 
isting around the new Dirac point energy. We show how- 
ever that, in triangular graphene superlattices (TGSs), 
there can be no states other than those of the new mass- 
less Dirac fermions around the energy of the new Dirac 
points. Therefore, doped or gated TGSs should provide 
a clear way to probe this new class of massless Dirac 
fermions that are absent in pristine graphene. 

A physical requirement for the discussed phenomenon 
is that the variation of the external periodic potential 
is much slower than the inter-carbon distance so that 
inter-valley scattering (between K and K') may be ne- 
glected 3 13 J and we limit our discussion to the low- 
energy electronic states of graphene which have wavevec- 
tors close to the K point. The Hamiltonian of the low- 
energy quasiparticles in pristine graphene in a pseudospin 



basis. 



Q j e*"^"" and y\j '^''^'^ (where ^ q j and ^ ^ 

are Bloch sums of 7r-orbitals with wavevector K on the 
sublattices A and B, respectivel y, a nd k is the wavevector 
from the K point), is given by 17 1 



Hq = hvQ {-ia^dx - icTydy) , 



(1) 



where vq is the group velocity and cr's are the Pauli ma- 
trices. The eigenstates and the energy eigenvalues are 
given by 



V2\ 



and 



E°{k) = shvok, 



(2) 



(3) 



respectively, where s = ±1 is the band index and 0k is the 
polar angle of the wavevector k. Equation ^ indicates 
that the pseudospin vector is parallel and anti-parallel 
to the wavevector k in the upper band (s = 1) and in 
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the lower band (s = —1), respectively. Moreover, the 
pseudospin vector is always parallel to the group velocity. 

Let us first consider the case that a ID potential V{x), 
periodic along the x direction with periodicity L, is ap- 
plied to graphene. The Hamiltonian H then reads 



H = hvQ {—icTxdx — io-ydy + I V{x)/?ivo) 



(4) 



where / is the 2x2 identity matrix. After a similarity 
transform, H' = u\HUi, using the unitary matrix 



where a{x) is given by (isj 



[a{x)/2 
»(x)/2 



a{x) 



2 / V{x')dx' /TiVQ, 

10 



we obtain [l£ 



H' = hvo 



-idx 

ia(x)Q 



,ia(x)Q 



(5) 



(6) 



(7) 



To obtain the eigenstates and energy eigenspectrum of 
H' in general, using a plane wave spinor basis set, we 
need an infinite number of plane waves with wavevec- 
tors different from one another by the reciprocal lattice 
vectors of the superlattice. (A reciprocal vector of the su- 
perlattice is given by G™ = m (27r/L) x = m Gq x where 
m is an integer.) However, if we are interested only in 
quasiparticle states whose wavevector k = p + Gm/2 is 
such that IpI <C Gq, we could treat the terms containing 
dy in Eq. d?]) as a perturbation since Gm is along x. H' 
may be reduced to a 2 x 2 matrix using the following two 
states as basis functions 



,»(p-G,„/2).r 



(8) 



[Note that the spinors ^ q ^ ^'^'^ 1 ) '^'^^ have a dif- 
ferent meaning from ^ q ^ 1 ) that were defined 

before because of the unitary transformation.] 

In order to calculate these matrix elements, we expand 

e^"(=") as 



(9) 



1 — — 00 



where /; \y] 's are coefficients determined by the periodic 
potential V{x). One important thing to note is that in 
general 



l/d<l, 



(10) 



which can directly be deduced from Eq. The physics 
simplifies when the external potential V{x) is an even 
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FIG. 1: Schematic diagram showing an equi-energy contour 
(ellipse) with E = hvoko + hvo m Go/2 of the newly generated 
massless Dirac fermions. The quasiparticle wavevector k, the 
generalized pseudospin vector (see text) c, and the group ve- 
locity vector Vg are represented by solid, dashed and dash- 
dotted arrows, respectively, for graphene in an even periodic 
potential. 



function. Then, /;[I^]'s are all real [20|. For states with 
wavevector k very close to Gm/2, the 2x2 matrix M 
whose elements are calculated from the Hamiltonian H' 
with the basis given by Eq. ^ can be written as 

M = hvo {px(Tz + fmPyCTy) + flVQ mGo/2 ■ I . (11) 

After performing yet another similarity transform M' = 
uImU2 with 



we obtain the final result: 



1 1 

-1 1 



(12) 



M' = hvQ {pxCTx + fmPyCTy) + flVQ m Go/2 ■ I . (13) 



The only difference of the Hamiltonian in Eq. from 
that in Eq. ([T]), other than a constant energy term, is 
that the group velocity of quasiparticles moving along 
the y direction has been changed from vq to fmVo 



211. 



Thus, the electronic states near k — Gm/2 are also those 
of massless Dirac fermions but having a group velocity 
varying anisotropically depending on the propagation di- 
rection. Moreover, the group velocity along the y direc- 
tion is always lower than vq [Eq. ((T0|) ] regardless of the 
form or magnitude of the periodic potential V{x). 

The eigenstate and the energy eigenvalue of the matrix 
M' are given by 



1 



V2 ( se'f^ 



(14) 



and 



Es{p) = shvo^pl + \fm?pl + mGo/2 , (15) 

respectively, where (p-p is the polar angle of the pseu- 
dospin vector c of 95s,p, which is parallel to s{pxX + 
fmPyy)- The spinor </3s,p, however, should not be con- 
fused with the one in Eq. ^ representing the sublattice 
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degree of freedom, or with the one in Eq. ([8]). A double 
prime in Eq. (jl4p emphasizes this point. 

The eigenstate il^s^ui^) of the original Hamiltonian H in 
Eq. (HI) can be obtained by using Eqs. dS]), dH), (HH) and 
()14|) . Since the unitary transforms conserve the inner- 
product between eigenstates, if a generalized pseudospin 
vector for the original Hamiltonian H in Eq. ^ is de- 
fined as the pseudospin vector of the transformed Hamil- 
tonian M' , i. e. , c, the scattering matrix elements be- 
tween states of these new massless Dirac fermions due to 
long-wavelength perturbations are described by the gen- 
eralized pseudospin in the same manner as those of the 
original massless Dirac fermions in pristine graphene are 
described with their pseudospin. 

On the other hand, the group velocity vector Vg is par- 
allel to s {PxX + f^PyV) [Eq. p5| ]. Therefore, in general, 
the three vectors p, c and are not collinear (Fig. [T]). 
However, it is obvious that if the wavevectors (p) of two 
electronic states are aligned or anti-aligned to each other, 
so are their generalized pseudospin vectors, as in pris- 
tine graphene, resulting in a maximum or a zero overlap 
between the two states, respectively. If V{x) is not an 
even function, the dispersion relation of the new mass- 
less Dirac fermions remains the same as Eq. ([T5|) but a 
generalized pseudospin vector may not be defined [23 |. 

Similarly, for graphene in slowly varying 2D periodic 
potential, new massless Dirac fermions are generated cen- 
tered around the wavevectors kc — G/2 where the G's 
are the superlattice reciprocal vectors. A state with 
wavevector k around kc mixes strongly with another 
state with wavevector k — G by the superlattice poten- 
tial. Applications of the same argument that we made 
use of in the case of ID graphene superlattices result in 
linear band dispersions. 

Even though new massless Dirac fermions are gener- 
ated in ID graphene superlattices, because there is no 
SBZ boundary perpendicular to the periodic direction, 
they are obscured by other states, and there is no new 
value of energy at which the density of states vanishes. 
In a 2D rectangular graphene superlattices, the SBZ is 
a rectangle. It turns out that the energy separation at 
the SBZ corners also vanishes due to the chiral nature 
of graphene 0. Therefore, in 2D rectangular graphene 
superlattices, again, there are states other than the new 
massless Dirac fermions in the range of the new Dirac 
point energy. However, as we show below, in TGSs, 
there can exist an energy window within which the only 
available states are the newly generated massless Dirac 
fermions. 

As an illustration, we consider a TGS shown in 
Fig. [2lja). The external potential is of a muffin-tin type 
with value Uq in a triangular array of disks of diameter 
d and zero outside of the disks. The spatial period of the 
superlattice is L. Figure [2I^b) shows the SBZ of a TGS. 

Figure OJc) shows the electron energy separation be- 
tween states in the first and the second band above the 




FIG. 2: (a); A TGS with muffin-tin type of periodic potential 
with a spatial period L. The potential is Uo inside the gray 
disks with diameter d and zero outside, (b) The SBZ of a 
TGS. (c) The energy separation AE between states in the first 
and the second band above the original Dirac point energy 
versus the wavevector k along the path KMK' in a TGS given 
by Uo =0.5 eV, L = 10 nm, and d = 5 nm. 

original Dirac point energy along the path KMK' in the 
SBZ [Fig.lHb)] for a TGS. The energy separation at the 
corner, or the K point, of the SBZ is largest, contrary 
to that of the rectangular graphene superlattices where 
the energy separation closes at the SBZ corners 0; but 
that at the M point is zero. New massless Dirac fermions 
are thus formed around the M points. With the set of 
potential parameters in Fig. [2] {Uq =0.5 cV, L = 10 nm, 
and d — b nm), the energy separation at the K point 
is 82 meV, much larger than room-temperature thermal 
energy. This energy separation can be tuned by changing 
the superlattice parameters. 

Figure [S^a) shows the energy dispersions of the first 
and the second band of the considered TGS. We can 
see the linear energy dispersion relation at the M points 
[Fig. m^c)]. Close to the original Dirac point energy 
(i? = 0), the density of states (DOS) varies linearly with 
energy, similar to that of pristine graphene, except that 
the slope is larger because of the reduced band velocity. 
At around = 0.16 eV, there exists another energy value 
where the DOS vanishes also linearly. 

In conclusion, we have shown that a new class of mass- 
less Dirac fermions are generated in graphene when a 
periodic potential is applied and we have studied the 
novel characteristics of these quasiparticles. Moreover, 
in triangular graphene superlattices, there can exist en- 
ergy windows where there are no other states than these 
new quasiparticles. The triangular graphene superlat- 
tices thus should provide a good platform for experimen- 
tal probing of the new massless Dirac fermions predicted 
here. 
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FIG. 3: (color online) (a): Energy dispersion relation of a 
TGS with external potential with Uo ~ 0.5 eV, L = 10 nm 
and d — 5 nm for the first and the second band above the 
original Dirac point energy as a function of wavevector k from 
the original Dirac point. Arrows indicate the M points of the 
SBZ around which new massless Dirac fermions are generated, 
(b): The DOS of charge carriers in electron orbits (bright and 
red) and hole orbits (dark and blue) in the TGS characterized 
in (a) . The original Dirac point energy is set at zero. Dashed 
black line shows the DOS of pristine graphene. The arrow 
indicates the new Dirac point energy. 
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Note added in proof. - After submission of this Letter, 
an angle-resolved photoemission experiment on graphene 
on Ir(lll) surface resulting in superlattice formation was 
reported [23j in which minigap openings at the SBZ 
boundary are found and evidence of replicas of the pri- 
mary Dirac cone observed. 
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M' 
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hvo 
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